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Abstract 

A parameter optimization framework has earlier 
been developed to solve the problem of partitioning a 
centralized controller into a decentralized, hierarchical 
structure suitable for integrated flight/propulsion control 
implementation. This paper presents results from the 
application of the controller partitioning optimization 
procedure to IFPC design for a Short Take-Off and 
Vertical Landing (STOVL) aircraft in transition flight. 
The controller partitioning problem and the parameter 
optimization algorithm are briefly described. Insight is 
provided into choosing various "user" selected parameters 
in the optimization cost function such that the resulting 
optimized subcontrollers will meet the characteristics of 
the centralized controller that are crucial to achieving the 
desired closed-loop performance and robustness, while 
maintaining the desired subcontroller structure constraints 
that are crucial for IFPC implementation. The 
optimization procedure is shown to improve upon the 
initial partitioned subcontrollers and lead to performance 
comparable to that achieved with the centralized 
controller. This application also provides insight into the 
issues that should be addressed at the centralized control 
design level in order to obtain implementable partitioned 
subcontrollers. 

Introduction 

Large interconnected systems often exhibit 
significant coupling between the various subsystems thus 
requiring an integrated approach to controller design. 
Short Take-Off and Vertical Landing (STOVL) aircraft are 
an example of such subsystems. In STOVL aircraft, the 
forces and moments generated by the propulsion system 
provide control and maneuvering capabilities for the 
aircraft at low speeds thus creating the need for Integrated 
Flight/Propulsion Control (IFPC) system design. An 
approach to IFPC design [1] is to design a centralized 
controller considering the integrated airframe and 
propulsion system with all its interconnections as the 
design plant. Although such an approach yields an 
"optimal" design since it accounts for all subsystem 
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interactions, it results in a high-order controller which is 
difficult to implement and validate. In aircraft design, it 
is the responsibility of the engine manufacturer to ensure 
that the propulsion system will provide the desired 
performance when installed in the airframe. The engine 
manufacturer therefore needs a separate engine controller 
to be able to independently perform extensive testing to 
assure adequate performance and integrity of the 
propulsion system in the presence of operational and 
safety limits. This requirement imposes the need for 
decentralized implementation of IFPC systems. One 
approach to direct decentralized design of IFPC systems 
is presented in Ref. [2], This approach consists of 
"partitioning" the overall system into loosely coupled 
subsystems and then performing a decentralized control 
design considering one subsystem at a time. Although 
such an approach results in low -order, independently 
implementable subsystem controllers (referred to as 
"subcontrollers"), it has the disadvantage that it does not 
easily account for all the interactions between various 
subsystems. 

An approach to IFPC design which combines the 
"best" aspects of the centralized and decentralized 
approaches is presented in Ref. [3]. This approach 
consists of first designing a centralized controller 
considering the airframe and propulsion systems as one 
integrated system, and then partitioning the centralized 
controller into decentralized airframe and engine 
subcontrollers with a specified interconnection structure. 
Here, "partitioning" means the process of approximating 
the high order centralized controller with two or more 
lower order subcontrollers, with a specified coupling 
structure, such that the closed -loop performance and 
robustness characteristics of the centralized controller are 
matched by the partitioned subcontrollers. A meaningful 
trade-off between subcontroller complexity and achievable 
performance for the integrated system can be performed 
by evaluating various controller partitionings of different 
levels of complexity against the performance baseline 
established with the centralized controller. 

The most suitable decentralized control structure 
for IFPC systems is hierarchical with the airframe (flight) 
controller generating commands for the aerodynamics 
control surfaces as well as for the propulsion subsystem. 
This hierarchical structure will be discussed in the next 
section. A stepwise approach to determining partitioned 



subcontrollers with the decentralized hierarchical structure 
from a centralized 1FPC design is presented in Ref. [4]. 
The procedure of Ref. [4], however, is ad-hoc in nature 
and can result in unacceptable degradation in closed-loop 
performance and robustness from that obtained with the 
centralized controller. A parameter optimization 
framework which can be effectively used to improve upon 
the subcontrollers obtained by the procedure of Ref. [4] so 
as to "more closely" match the centralized controller 
closed-loop performance ami robustness characteristics 
was discussed in Ref. [5]. The application of the 
partitioned subcontroller parameter optimization algorithm 
to IFPC design for a Short Take-Off and Landing (STOL) 
aircraft was also presented in Ref. [5]. The results 
obtained for the STOL aircraft application demonstrated 
the potential of the parameter optimization approach to 
obtain partitioned subcontrollers that match the 
performance and robustness characteristics of the 
centralized controller while maintaining the desired 
subcontroller structure. The objective of this paper is to 
present the results from application of the controller 
partitioning parameter optimization procedure to IFPC 
design for a STOVL aircraft. The STOVL aircraft 
problem presents a significant increment in complexity 
over the STOL example studied earlier. Extensive insight 
is to be gained about the challenges associated with 
partitioning a complex IFPC centralized controller by 
application to the STOVL example. 

The paper is organized as follows. Hie controller 
partitioning problem and the optimization framework is 
first briefly described. The STOVL vehicle model and the 
centralized as well as partitioned IFPC controller structure 
is then described along with a brief discussion of the 
performance results obtained with the initial partitioned 
controllers derived by using the procedure described in 
Ref. [4]. The results from application of the parameter 
optimization procedure are then presented and 
performance comparisons are made between the 
centralized controller, initial partitioned subcontrollers and 
the optimized subcontrollers. Issues related to the 
structure of the propulsion subsystem controller for IFPC 
implementation are then discussed in light of the results 
obtained by the parameter optimization procedure. 

Controller Partitioning: 

Problem and Optimization Algorithm 

The decentralized, hierarchical controller 
partitioning structure is shown in Fig. 1 where the 
subscripts and superscripts "a" and "e" refer to airframe 
and propulsion system (engine) quantities, respectively, 
subscript "c" refers to commands, and the variables z are 
the controlled outputs of interest with e being the 
corresponding errors. The intermediate variables^ 
represent propulsion system quantities that affect the 
airframe, for example propulsion system generated forces 
and moments. 

The controller partitioning problem of Fig. 1 can 
be stated as follows: 


Given: A centralized controller K(s) s.t. 


Find: 


u(s) = K(s) 


e(s) 

|y(s). 
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and a particular set of the interface variables 


Decentralized airframe and engine 
subcontrollers, K*(s) and K e (s), respectively, 
with 


u,(s) 

s {s) . 

So that: 


= K*(s) 


e,(s) 

y.(s) 


, and u e (s) = K'(s) 


e„(s) | 
e e (s) 

y c (s) 


( 2 ) 


The closed-loop performance and robustness 
with the subcontrollers K B (s) and K c (s) match 
those with the centralized controller K(s) to a 
desired accuracy. Furthermore, the engine 
subcontroller K e (s) should have the structure of 
a command tracking controller for the interface 


variables z^ to allow for independent check- 
out of the propulsion system. 

The approach for solving the controller partitioning 
problem, as discussed in Ref. [5], consists of optimization 
of a suitable cost function over the state-space parameters 
of the partitioned subcontrollers using an analytical 
expression for the gradient of the cost function. The 
initial partitioned subcontrollers to start the parameter 
optimization are obtained using the procedure discussed in 
Ref. [4] and the cost function to be minimized is chosen 
to be of the form 


J (p) = yp) + j Mck (p> 

where p is the vector of parameters over which the 
optimization takes place, J^p) reflects the performance 
requirements (including robustness) and J tack (p) reflects 
the z^ tracking requirement. The performance cost is 
chosen to be 


3 tJv) = fe (o k [W i (jw)-(K(ju3)-K(p)(jio)) -W^uo)]) 2 du: 
®, k 

(3) 

where o k [*] denotes the k" 1 singular value of a matrix, and 

K(s) is the "equivalent" centralized controller obtained by 
assembling the partitioned subcontrollers. Plant 
information consisting of the state -space representation of 

the transfer function matrices G(s) and 6^(s), defined by 
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V. 




= G(s)u ; z =<j_(s)u 


(4) 


is used to obtain the "equivalent" centralized controller, 


K(s), as described in Ref. [5]. The above choice of 
Jpeii^P) corresponds to the H 2 norm of the weighted 
difference between the designed centralized controller and 
the "equivalent" centralized controller. The frequency 
band, is the interval over which a good match 

between the two controller transfer function matrices is 
sought while the input and output frequency weightings, 
W,(s) and W 0 (s), respectively, allow for emphasizing 
certain frequency ranges and directions in obtaining this 
good approximation. 

The z^ tracking cost J tflci: (p) is chosen to be 

= IE V(rrL(i«»-t '(pXW d«*> (5) 

i ■' 

where: T 1 (s) is the transfer function vector from the 
airframe commands z t to the i ,h interface variable zj with 

the centralized controller, T *(s) is the transfer function 
vector from the airframe commands z t to the I th 

commanded interface variable^ 1 , with the partitioned 
airframe controller; \ is a weighting which could be 
chosen to be frequency dependent; and H 2 denotes the 
Euclidean norm of a row vector. Note that for the 
partitioned subcontrollers to closely match the performance 
achieved with the centralized controller, the response of 
the interface variables z w to airframe commands z a with 
the partitioned subcontrollers must match the 
corresponding response with the centralized controller 
because the interface variables (e.g., propulsion system 
generated thrusts) significantly affect the airframe 

responses z a . Thus requiring T ^(s) to closely match 

f l (s), as reflected in J lclck (p), will result in partitioned 
controllers such that z ^ appears to be a command for 

z The weights Xj allow the control designer to provide 
relative weighting for enforcing the command tracking 
structure among the various elements of z ^ and for 
trading off the performance cost against the cost of 
enforcing this command tracking structure. 

In Ref. [5], a procedure was presented for developing 
the analytical gradients of the cost functions described in 
(3) and (5) with the parameter vector p consisting of the 
matrix elements in the state-space realization of the 
subcontroller transfer function matrices K*(s) and K c (s). 
Based on these analytical expressions for the gradients, a 
numerical algorithm for minimizing the cost J(p) was 


developed. This algorithm is shown in flowchart form in 
Fig. 2 and is briefly described in the following. 

The fixed data used by the algorithm of Ref. [5] are 
state-space representations for the plant system matrices 

G(s) and (^(s) (as defined in (4)), the centralized 
controller K(s), the weighting matrices W,(s) and W o (s), 
the weights X^ and the partitioning structure consisting of 
subcontroller inputs and outputs and the interface 
variables. The initial partitioning is obtained using the 
stepwise procedure described in Ref. [4] or as a result of 
order reduction on an earlier optimized partitioning. This 
initial partitioning is converted to a "minimal parameter" 
form described in Ref. [5] and an initial parameter vector 
is generated. For a given parameter vector, the cost and 
the analytical expression based cost gradients are then 
calculated. The Broyden-Fletcher-Goldfarb-Shanno 
method (see Ref. [6]) for search direction and linesearch 
is used to update the parameter vector such that both the 
cost function and the norm of the gradient vector is 
reduced. This linesearch is constrained so as to maintain 
the stability of the subcontrollers. Various convergence 
criteria involving reduction in the cost function between 
successive parameter updates, the norm of the difference 
between previous and updated parameter vector, and the 
norm of the gradient vector are used to check for 
"optimality". The output of this procedure is a set of 
partitioned subcontrollers which can be analyzed for 
performance and stability robustness. If these are 
unacceptable, then the optimization can be continued with 
either tighter convergence criteria or modified choice of 
the weighting matrices W,(s) and W 0 (s) and/or the weights 
X, 

STOVL Aircraft Model and 
Initial Controller Partitioning 
The controller partitioning optimization procedure 
discussed above was applied to the IFPC design for a 
STOVL aircraft in the decelerating transition during 
approach to hover landing flight phase. A schematic 
diagram of the aircraft is shown in Fig. 3. The aircraft is 
powered by a two-spool turbo fan engine and is equipped 
with the following control effectors: left and right elevons 
used collectively as elevator and differentially as ailerons; 
rudder, ejectors to provide propulsive lift at low speeds 
and hover, a 2D-CD (two dimensional convergent- 
divergent) vectoring aft nozzle; a vectoring ventral nozzle 
for pitch control and lift augmentation during transition; 
and jet reaction control systems (RCS) for pitch, roll and 
yaw control during transition and hover. Engine 
compressor bleed flow (WB3) is used for the RCS 
thrusters and the mixed engine flow is used as the primary 
ejector flow. The aircraft and engine model and the 
design of the centralized controller for a linear integrated 
design model are discussed in detail in Refs. [7,8]. The 
centralized controller was partitioned into decoupled lateral 
and longitudinal-plus-engine subcontrollers as discussed in 
Ref. [8], and the longitudinal-plus-engine controller was 
further partitioned into separate airframe and engine 
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subcontrollers which have the decentralized, hierarchical 
structure of Fig. 1. In the following, the vehicle model is 
first summarized, and the performance results with the 
initial partitioned subcontrollers obtained from Ref. [4] are 
briefly discussed. 

The linear integrated aircraft longitudinal dynamics and 
engine dynamics small perturbation model is of the form 

x = Ax + Bu (6) 

where the state vector is 

x = [N2,N25,Tmhpc,Tmpc,Tmhpt,Tmlpt,u,w,q,04i] T 
with 

N2 * Engine Fan Speed, rpm 

N25 = High Pressure Compressor Speed, rpm 

Tmhpc = High Press. Compressor Metal Temp., °R 

Tmpc - Burner Metal Temp., °R 

Tmhpt = High Pressure Turbine Metal Temp., °R 

Tmlpt = Low Pressure Turbine Metal Temp., °R. 

u = Axial Velocity, ft/s 

w = Vertical Velocity, ft/s 

q = Pitch Rate, rad/s 

0 = Pitch Attitude, rad 

h = Altitude, ft 

The control inputs partitioned into airframe and engine 
control inputs are 

H. = [6e,AQR,ANG79,ANG8] T 
5 = [WF,A8,ETA,A78 ] t 

with 

be = Elevator Deflection, deg 

AQR = Pitch RCS Area, in 2 

ANG79 - Ventral Nozzle Vectoring Angle, deg 

ANG8 * Aft Nozzle Vectoring Angle, deg 

WF = Fuel Flow Rate, lbm/hr 

A8 = Aft Nozzle Area, in 2 

ETA * Ejector Butterfly Angle, deg 

A78 = Ventral Nozzle Area, in 2 

The controlled outputs for the airframe and engine systems 
are 

Z. = [Vv,Qv,y] T ; z. = N2 

where Vv=V+0.1V, Qv=q+O.30 with 

V = True Airspeed, ft/s 

V = Acceleration Along Flight Path, ft/s 2 

y = Flight Path Angle, deg 

and the other outputs as discussed under state description 
with units of q and 0 in degrees. As discussed in Ref. [7], 
the above choice of z n corresponds to providing the pilot 
with an acceleration command, velocity hold system in the 
forward axis; pitch rate command, attitude hold system in 
the pitch axis; and direct command of the flight path angle 
for vertical axis control. The choice of z e allows for 
setting the engine operating point independent of the 
aircraft maneuver. 


The inputs to the airframe and the engine controllers are 
the tracking errors e a and e c corresponding toz^ and z e 
respectively, and the measurement feedbacks 

y, - [V,V,0,q,y] T ; y c = [N2,WB3] T 
where WB3 is the compressor bleed flow demanded by 
the RCS control. 

The interface from the propulsion system model to the 
airframe model is defined by the gross thrusts from the 
three engine nozzle systems, i.e. 

z m = [FG9,FGE,FGV] T 

where 

FG9 * Aft Nozzle Gross Thrust, lbf 
FGE = Ejector Gross Thrust, lbf 
FGV = Ventral Nozzle Gross Thrust, lbf. 

Using the procedure described in Ref. [4], an initial 
partitioning of the centralized longitudinal-plus-engine 
controller into separate airframe and engine subcontrollers 
was obtained. These initial partitioned subcontrollers 
provided close matching of the partitioned closed-loop 
system response to the centralized closed -loop system 
response for aircraft velocity and flight path commands 
(Vv c and y c ) and the engine fan speed command (N2 C ). 
However, this partitioned system response to the pitch 
variable command (Qv c ) showed significant deviation from 
the corresponding response for the centralized system. 
Shown in Fig. 4 is the closed-loop response of the system 
with the centralized and initial partitioned subcontrollers 
for an example pitch variable command. Although the 
level of decoupling in the Vv, y and N2 response to Qv c 
with the initial partitioned subcontrollers is comparable to 
that obtained with the centralized controller, there is 
significant degradation in the tracking of the Qv command 
itself. The partitioning parameter optimization algorithm 
was applied to this example to investigate whether the 
response to pitch variable command can be improved. 

Controller Partitioning Optimization 
As mentioned earlier, the controller partitioning 
parameter optimization algorithm was earlier applied to a 
linear model of a STOL (Short Take-Off and Landing) 
aircraft and the results were presented in Ref. [5]. The 
STOVL problem being addressed in this paper is 
considerably more complex than the STOL example in 
that the centralized controller is not square (K(s) has 
dimension 8 by 10 for the STOVL problem as opposed to 
4 by 4 for the STOL example), and there are 3 interface 
variables (z^) as opposed to just one for the STOL 
example. Further insight is to be gained into the 
suitability of the controller partitioning parameter 
optimization procedure and the proper selection of the 
various weighting factors by applying the procedure to this 
complex STOVL example. In the following, the issues 
related to the proper choice of weighting factors are first 
discussed and then the results are presented for a set of 
optimized subcontrollers. 
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Choice of Weighting Factors in the Optimization Cost 
Initially, the weighting matrices for the performance cost 
function J^p) were chosen to be consistent with those 
that led to the successful results for the STOL example, 
i.e., W,«I, and W H (s) = G(s)-[I+K(s)G(s)]-\ where I is an 
appropriately dimensioned identity matrix and G(s) is the 
plant system matrix defined earlier. This choice of 
weighting then corresponds to minimizing the H 2 norm of 

(K(ju>)-K(p)(j<*>)) G(jw)-[I+K(ju))G(ju))]'' (7) 

which is the frequency weighted loop transfer matrix error 
at the control inputs u . The weights \ in J^lp) were 
chosen to be 

= tirLd^UU-' (8) 

which corresponds to normalization of the tracking error 
for I th interface variable by the euclidean norm of the 
steady-state response of the i 01 interface variable to the 
airframe commands (z a ) for the centralized system. 

Exercising the optimization algorithm with the above 
choice of weightings resulted in significant reduction of 
the total cost, J(p), as well as reduction in both the 
individual elements of the cost, J^p) J tnC k(P)- 
However, when the closed-loop system was analyzed with 
the optimized partitioned subcontrollers, the performance 
was much degraded over the initial partitioning. Note that 
the choice of W D (s) corresponding to (7) was driven by the 
small gain theorem according to which the closed-loop 
system with the partitioned subcontrollers will remain 

stable if |(K0^)-K(p)0to))G0^) , [I + K(jto)G(j^)]* 1 L < 1> 
where denotes the EL norm. This type of weighting 
is used in controller order reduction problems (see Ref. 
[9]) and was successful for the STOL controller 
partitioning example because the controller and the plant 
were square. For a square system, matching the loop 
transfer function matrix at the control inputs (u) will in 
general imply a good match with the loop transfer function 
matrix at the controlled outputs ( z ), which corresponds to 
matching closed -loop performance. This sort of 
relationship does not necessarily hold true for a non- 
square system as was experienced for the STOVL 
example. 

Next, the weights in the performance cost, J^p) , were 
modified such that Wj(s)=G z (s) and W 0 (s)=l, where G z (s) 
is defined such that z=G z (s)u. This choice of weighting 
corresponds to minimizing the H 2 norm of 

G z (jio) <K(jto) -K(jto)> (9) 

which is the loop transfer matrix error at the controlled 
outputs z. The "optimal" partitioned subcontrollers 
obtained by exercising the optimization algorithm with 
these modified weightings were such that these led to 
improved closed-loop performance over the initial 
partitioned subcontrollers in terms of more closely 
matching the decoupled command tracking properties of 
the system with the centralized controller. However, these 


optimized subcontrollers had excessive control and control 
rate requirements which will result in significant system 
performance degradation when implementing the 
subcontrollers with control actuation and rate limits. 
Furthermore, the required interface variable (z^) 
command tracking structure for the engine subcontroller 
was not retained by these optimized subcontrollers. 
Detailed analyses of the above results revealed that the 

equivalent controller JC(s) obtained by assembling these 
optimized subcontrollers has totally different input/output 
response characteristics as compared to the designed 
centralized controller K(s). Note that for the STOVL 
example there are only four controlled variables, z, 
whereas there are eight control inputs, u . Due to this 
redundancy in control effectors, theoretically there are an 
infinite number of solutions for a control law which 
matches the loop transfer function matrix at the controlled 
outputs with the designed centralized controller. Without 
any constraints on the control usage, "larger" control 
provides "better" matching of the loop transfer function 
matrix, and so the parameter optimization algorithm tends 

to a solution for K(s) which has excessive control 
requirements. Further note that the formulation for the 
"tracking” cost, J tI>ck (p)> was based on the hypothesis that 
"for the partitioned subcontrollers to closely match the 
performance achieved with the centralized controller, the 
response of the interface variables to airframe 
commands z a with the partitioned subcontrollers must 
match the corresponding response with the centralized 
controller". This hypothesis is no longer true, for the 
above choice of weightings in due to the controller 
redundancy effect discussed earlier. Although the choice 
of J twck forces the z^ response to z t command input with 

the optimized partitioned subcontrollers to match thez^ 
response with the centralized controller, the z w response 
with the partitioned subcontrollers is itself very different 
from the corresponding response with the centralized 
controller. Thus, for the optimized partitioned 
subcontrollers the z^ and z^ responses are very different 
from each other and the desired command tracking 
relationship between z a and z a is not maintained. 

In order to overcome the problems due to control 
redundancy and to force the "same" control solution for 
the optimized subcontrollers as for the centralized 
controller, the optimization was performed with Wj(s) and 
W 0 (s) chosen to be appropriately dimensioned identity 
matrices, i.e. minimize the H 2 norm of the unweighted 
controller approximation error 

(K(jtn) ”K(ju))) . (10) 

With this choice of the optimized subcontrollers 
showed improvements in all areas in which there were 
problems with the previous formulations of J^, however 
the resulting engine subcontroller was very high bandwidth 
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with virtually flat response across a large 

frequency band for the closed-loop engine subsystem. The 
parameter vector over which the optimization was 
performed included elements of the direct feedthrough 
matrix from the error in interface variable command 
tracking (e^ -z^ -z tt ) to the engine control inputs (u e ) for 
the engine subcontroller. As the optimization proceeded, 
these elements tended to become large thus implying a 
large bandwidth for the interface variable command 
tracking portion of the engine subcontroller. In order to 
limit the engine subcontroller bandwidth and provide 
adequate control gain attenuation, the engine subcontroller 
structure was modified to keep this direct feedthrough 
matrix zero for the partitioning optimization. To provide 
further frequency roll-off, the weights \ in were 
modified to be 

(ii > 

The results with this choice of weighting factors in 
and J lnd: are discussed in the following subsection. 

Optimization Results and Discussion 

Shown in Fig. 5 are the maximum singular values (5) 
of three controller partitioning error measures 
corresponding to (10), (9) and (7) for the initial and the 
optimized partitioned subcontrollers. The optimized 
partitioning shows improvement over the initial 
partitioning for all three measures although 
corresponded to only minimizing H 2 norm of (10). When 
the optimization was done with corresponding to 
either (7) or (9), the other two error measures ((10) and 
(9) or (7), respectively) showed an increase over a 
significant frequency region leading to the problems 
discussed earlier. The significant reduction in the 
unweighted controller approximation error 

(a[K(Jto)-K(jco)]) shows that the optimized partitioned 
subcontrollers better match the centralized controller, the 
significantly decreased error in the loop transfer function 

matrix at the controlled outputs ( o[G z 0^MK(ju)) -K(jio))] ) 
indicates that the optimized subcontrollers will better 
match the decoupled command tracking properties of the 
centralized controller, and the reduced error in the 
weighted loop transfer function matrix at the control inputs 

(5[(K.(Ju3)-K(ju)))<j(jco)[l+K(jco)G(jtJO)]" 1 ]) indicates that 
the optimized subcontrollers will better match the stability 
robustness characteristics of the centralized controller with 
respect to uncertainties at the control inputs. Note that for 

the optimized subcontrollers, |(K-K)G[I+KG]~ , L < 1 , 
which guarantees that the closed-loop system will be 
stable with the optimized partitioned subcontrollers. 

The closed-loop response to pitch variable command for 
the optimized partitioned subcontrollers is shown in Fig. 
4. The optimized partitioned subcontrollers provide the 
desired improvement over the initial partitioning in the Qv 
command tracking response while the level of decoupling 
in the Vv, y and N2 responses is maintained. For all the 


other command inputs (z c ) also, the optimized 
subcontrollers provided equally good or slightly improved 
performance as compared to the initial partitioned 
subcontrollers. The control input (u) requirements for all 
commands with the optimized subcontrollers were quite 
similar to those with the centralized controller. 

The results in Fig. 6 show the effect of J lnick in imposing 
the engine subcontroller command tracking structure 
discussed earlier. Shown in Fig. 6 are the frequency 


response magnitude plots for the three z ^ — z ^ 
responses with the initial and optimized engine 
subcontroller. The optimization results in increased 
tracking bandwidth for aft nozzle and ejector thrust 
commands (FG9 C and FGE C ), but a decrease in the 
tracking bandwidth for ventral nozzle command (FGV C ) 
accompanied with a somewhat increased steady-state 
tracking error. Also, although not shown here, there was 
significant coupling in the FGV response from the FG9 
and FGE commands indicating that this optimized 
subcontroller will not meet the requirements for 
independent check-out of the propulsion subsystem. This 
shortcoming of the optimized engine subcontroller can be 
overcome by varying the weightings Aj in J tnck as was 
successfully demonstrated in Ref. [5] for the STOL 
example. However, this was not done for the current 
STOVL example because experience gained by exercising 
the partitioning optimization algorithm on this complex 
problem suggests a modification in the problem 
formulation which will be discussed next. 

Consider the engine subcontroller as consisting of two 
components: 


u(s) = K * (s)e (s) + K e ($i 


e c (s) 

y c (s) 


( 12 ) 


The centralized controller, K(s), contains a sub-block 


corresponding to the same inputs and outputs as K c e (s), 
the z e command tracking and feedback augmentation 
portion of the engine subcontroller. The choice of 
directly provides a "design" constraint on K e e (s) in terms 
of closely matching the equivalent portion of the 
centralized subcontroller. There is no such direct 
information available from the centralized controller 


regarding the design of K ^(s), the z^ command tracking 
portion of the engine subcontroller. The centralized 
controller only provides guidelines on the minimumz^ 
command tracking bandwidth required of the engine 
subcontroller. As discussed in Ref. [4], there are other 
constraints placed on the design of K ^(s) such as 
disturbance rejection requirements, robustness to modelling 
uncertainties and control actuation limits etc. For the 


initial controller partitioning [4], a K^s) is designed 
such that it best meets the various design requirements. 
However, the K ^(s) gets modified during the partitioning 
optimization, and since the J inck formulation does not 
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adequately reflect all the design requirements for K ^(s), 
the optimization can result in an unacceptable K '^(s). An 
approach to keeping the K c ea (s) portion of the engine 
subcontroller fixed during the optimization is currently 
being investigated. 

In implementation of propulsion system control laws, 
the command for the engine operating point is determined 
from an open-loop schedule which is based on some 
measure of the gross thrust demanded from the propulsion 
system. For the STOVL example, this corresponds to the 
fan speed command being a function of the gross thrust 
commands for the three nozzles, i.e. N2 C = 
f(FG9 c ,FGE c ,FGV c ). Generalizing, this implies that the 
engine commands z e might be a function of the interface 
variable commands z a . In order for the engine 
subsystem to track the interface variable commands in the 
presence of such an outer loop, it is necessary that the 
engine subcontroller be such that it provides decoupling of 
the interface variable response to the engine commands, 
i.e. z c -*z^ * 0. Shown in Fig 7 is the response of the 

three interface variables, FG9, FGE and FGV, to a step 
engine command (N2 C = 200 rpm) for the centralized 
controller, the initial partitioned subcontrollers and the 
optimized partitioned subcontrollers. All the quantities 
shown in Fig. 7 correspond to perturbations from a trim 
condition. The initial partitioned engine subcontroller was 
designed to take into account the z e — » ► * 0 requirement 

aid this is reflected in the responses shown in Fig. 7. Thez c -+ z a 
decoupling requirement was not imposed in the centralized 
control design because it was thought that since the gross 
thrusts (zj significantly affect the airframe response (z ft ), 
imposing the criterion of decoupled command tracking for 
airframe and engine commands (z a and z e ) will result in 
the decoupling of thrust response to engine commands. 
However, due to the effect of control redundancy 
discussed earlier, the centralized controller is such that it 
results in significant coupling from the engine commands 
to the thrust response while decoupling of the airframe 
response is maintained by appropriate usage of the 
airframe control inputs (u a ). Since the optimized 
partitioned subcontrollers match the control usage of the 
centralized controller, via the formulation discussed 
earlier, these too result in significant coupling of the thrust 
response to fan speed commands as seen from Fig. 7. 

An approach to impose this z e — ► decoupling 
requirement in the parameter optimization procedure is to 
add an appropriately formulated cost function to the total 
cost J to be minimized. This approach is currently being 
investigated. Another approach might be to directly 
consider this decoupling requirement at the centralized 
control design level. However, this latter approach will 
require further research in control theory as most 
multivariable control design techniques such as Linear 
Quadratic Gaussian, H-, etc. do not have any direct means 


of incorporating design criterion which corresponds to 
penalizing an individual input/output response. 

Conclusions 

Results were presented from the application of a 
controller partitioning parameter optimization algorithm to 
Integrated FUght/Propulsion Control (IFPC) design for a 
Short Take-Off and Vertical Landing (STOVL) aircraft. 
Insight was provided into the effect of various user 
selected weighting parameters in the optimization cost and 
it was shown that with an appropriate choice of these 
weighting parameters, partitioned subcontrollers could be 
obtained that closely matched the closed -loop performance 
and robustness characteristics of the centralized controller. 
However, the current optimization problem formulation 
was found to be inadequate in terms of meeting the 
requirements placed by the need to be able to 
independently check-out the propulsion subsystem. This 
requirement is primarily that the engine subcontroller have 
a decoupled command tracking structure for the interface 
variables from the engine to the airframe, i.e. engine 
developed thrusts. A modification to the optimization 
procedure, which consists of keeping a subportion of the 
engine subcontroller fixed during optimization, is currently 
being investigated. As a result of this application study, 
a need was identified for developing modifications to 
modem multivariable control design techniques which will 
allow for penalizing an individual closed-loop input/output 
response in a multivariable system. 
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Centralized Control Loop 



Decentralized, Hierarchically partitioned control loop 


Fig. 1 Controller partitioning structure 



Fig. 4 Closed-loop system response of controlled 
variables (z) to Qv c 



Fig. 2 Flowchart for partitioning optimization algorithm 
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Fig. 3 Schematic of E-7D aircraft showing various control 
effectors 
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Fig. 5 Maximum singular value of various error measures 
for partitioned subcontrollers 
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